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I.  INTRODUCTION 


As  a  major  part  of  its  impact  acceleration  research  program, 
the  Naval  Biodynamics  Laboratory  is  investigating  the  dynamic  response 
of  the  head/neck  system  as  a  function  of  motion  and  anthropometric 
parameters.  One  goal  of  this  program  is  the  development  of  an  impact 
acceleration  injury  prediction  model  which  permits  reliable  inferences 
about  injury  probability  as  a  function  of  dynamic  and  physical  variables. 
Two  previous  Desmatics  reports  [2,3]  have  discussed  procedures  for 
incorporating  various  sources  of  auxiliary  information  into  a  model. 

This  report  describes  a  Monte  Carlo  investigation  of  the 
procedure  for  using  preinjury  "side  effect"  data  based  on  evoked 
potential  response.  The  notation  and  discussion  in  this  report  is 
a  continuation  of  that  of  the  previous  reports.  Thus,  the  reader 
may  find  a  review  of  those  reports  helpful  as  a  background. 


II.  THE  MONTE  CARLO  STUDY 


i 

i 

In  order  to  assess  the  performance  of  the  procedure  for  relatively 
small  samples,  500  Monte  Carlo  simulations  were  performed  using  a  sample 
size  of  25  observations.  The  goal  of  this  Monte  Carlo  investigation  was 
to  obtain  an  indication  of  the  effect  of  the  size  of  the  correlation 
coefficient,  p,  between  the  injury  tolerance  and  the  preinjury  side  effects 
on  the  mean  square  error  (MSE)  of  the  weighted  least  squares  estimates 
of  0.  Furthermore,  this  Monte  Carlo  study  assumed  no  prior  information 
linking  the  parameters  of  and  0^  *n  order  to  impose  a  more  stringent 
test  of  the  contribution  of  the  empirical  auxiliary  Information  to  de¬ 
creasing  the  MSE's  of  the  probit  estimates  of  the  elements  in  0^  when 
the  number  of  observations  is  small. 


A.  PROCEDURE 

Root  mean  square  error  (RMSE),  the  square  root  of  MSE,  was  used  as 
the  basis  for  performance  evaluation.  RMSE's  of  two  different  weighted 
least  squares  estimates  of  J3^  were  compared  for  five  values  of  p:  -0.1, 
-0.3,  -0.5,  -0.7,  -0.9.  The  first  estimate  of  £  was  derived  from  a  probit 
model  that  is  not  conditional  on  the  linear  regression  residuals.  This  is 
referred  to  as  the  "standard"  model.  The  second  estimate  of  was  derived 
from  a  probit  model  that  jl£  conditional  on  the  linear  regression  residuals. 
This  model  is  referred  to  as  the  "modified"  model. 


The  "standard"  and  "modified"  probit  models  and  the  linear  regression 
model  all  have  the  same  predictor  variable.  Values  for  the  predictor 
variable  are  sampled  from  a  uniform  distribution  between  -1  and  1.  For  each 
simulation  trial,  a  prior  unbiased  estimate  of-^^  was  sampled  from  a  normal 
distribution  with  mean  0^  and  standard  deviation  0.5.  This  estimate  of 
was  incorporated  into  both  the  "standard"  and  "modified"  problt  models 
as  outlined  in  [2]. 

Thus,  the  Information  employing  the  "standard"  model  can  be  represented 
by 

1  ■  6u  +  v 

yu  •  »(8oi  +  +  eli>  i  •  l . 25 

where  E(v)  *  0,  Var(v)  »  0.5, 

E(eli)  -  0,  Var<£11)  -  p^l-pj),  p£  -  *(801  +  B^j), 

and  r  is  an  unbiased  estimate  of  B^. 

Similarly,  the  information  conditional  on  the  linear  regression  residuals  can 
t-  represented  using  the  "modified"  probit  model  by 

yli  *  (l-p2)-1*  (B01  +  BjjXi  -pSt  )]  +  ,  i-1,  ...,  25 

where  EC^)  -  0,  VarC^)  -  Pjd-pp,  Pi  -  *[ (l-p2)_ls  (BQ1  +  Bj^  -  pSj)] 
th 

and  is  the  1  standardized,  maximum  likelihood,  linear  regression  residual. 

The  parameters  of  the  linear  regression  model  used  in  the  simulation 
were  801  »  B-,,  ■  O-  •  1.0.  The  starting  values  for  the  probit  iterations 


r 

for  8  ,  and  p  were  0.5,  0.5,  and  -0.5,  respectively. 

B.  RESULTS 

Figure  1  contains  the  Monte  Carlo  estimates  of  the  means,  standard 
deviations,  RMSE's  and  asymptotic  standard  deviations  for  the  estimates 
of  8^  and  8^  for  both  the  "standard"  and  "modified"  models.  Figure  1 
also  lists  the  Monte  Carlo  estimates  of  the  ratio  of  the  corresponding 
RMSE's  of  the  "standard"  model  to  the  "modified"  model  for  both  8q^  and 
8^.  Figure  2  lists  the  Monte  Carlo  estimates  of  the  mean,  standard 
deviation,  RMSE  and  asymptotic  standard  deviation  for  each  estimate  of  p. 
For  labeling  purposes,  estimates  of  Bq^  and  8^  from  the  "standard" 

A  A 

model  will  be  denoted  by  8^  and  8^,  whereas  estimates  of  8q^  and  8^^ 

A  ft 

from  the  "modified"  model  will  be  denoted  by  8^  and  B-q.  As  can  be 

ft  ft 

seen  from  Figure  1,  Bg^  and  8^  each  have  smaller  estimated  RMSE's  than 

A  A  ^  A 

their  counterparts  Bq^  and  B^.  (However,  the  bias  of  8qi  and  8^  is 

A  A 

slightly  larger  than  the  bias  for  Bq^  and  8^»  respectively.)  Increases 

in  the  absolute  value  of  p  tend  to  reduce  the  corresponding  RMSE's  of 

A  *  A 

B01  and  8^^,  with  the  tendency  more  pronounced  for  8jj. 

A 

Except  for  the  estimate  of  8^  for  p  ■  -0.7,  the  average  asymptotic 
standard  deviations  are  smaller  than  the  corresponding  average  estimated 
standard  deviations  derived  from  the  500  Monte  Carlo  trials.  This  is  true 
for  both  the  "standard"  and  "modified"  models. 

Figure  2  shows  the  Monte  Carlo  results  for  the  estimation  of  p  via 
the  "modified"  model.  The  average  estimated  p  has  a  bias  of  roughly  a 
little  over  0.1.  It  is  therefore  suggested  that  the  approximate  unbiased 
estimate  of  p  as  discussed  in  [3]  also  be  computed  in  addition  to  the 
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Parameter  Estimate  Standard  Deviation 


Parameter 
p  »  -0.1 

p  -  -0.3 
p  -  -0.5 
p  *  -0.7 
p  -  -0.9 


Figure  2: 


Parameter 


Estimate 

Standard  Deviation 

0.0062 

0.2917 

-0.1864 

0.2754 

-0.3741 

0.2422 

-0.5757 

0.1836 

-0.7554 

0.1333 

RMSE 

Asymptotic 
Standard  Deviation 

0.3064 

0.2497 

0.2979 

0.2429 

0.2730 

0.2213 

0.2217 

0.1820 

0.1967 

0.1319 

Monte  Carlo  Estimates  for  p  (correlation  coefficient)  Based 
on  500  Simulations. 
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estimate  of  p  from  the  "modified"  model.  Note  that  the  RMSE  of  the 
estimate  of  p  tends  to  decrease  with  increasing  absolute  values  of  p 
and  that  the  estimated  asymptotic  standard  deviation  is  always  smaller 
than  the  estimated  standard  deviation. 

C.  SUMMARY 

Auxiliary  empirical  information  based  on  a  preinjury  "side  effect" 
measurement  should  be  helpful  in  reducing  the  MSE  of  the  probit 
prediction  model  parameter  estimates  in  the  estimation  situation 
presented  in  the  previous  section.  Even  though  the  sample  size  was  only 
25  and  there  was  no  prior  information  linking  8  and  j^,  the  auxiliary 
empirical  information  noticeably  reduced  the  MSE  of  the  probit  estimate 
of  .  Had  there  been  prior  information  linking  jj  and  j^,  e.g., 
f^ll  *  ®12^°2’  a  muck  larger  reduction  in  MSE  would  have  taken  place. 
Fortunately,  this  parameter  constraint  situation  as  discussed  in  [3] 
seems  reasonable.  Therefore,  further  consideration  of  the  existence  of 
such  a  constraint  should  be  taken  into  account. 
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III.  DISCUSSION 


An  Important  statistical  topic,  somewhat  connected  with  the 
simultaneous  analysis  of  the  injury  and  preinjury  data,  is  the  subject 
of  "seemingly  unrelated"  regression  estimation.  (See  [1],  for  example.) 

By  "seemingly  unrelated"  regression  estimation,  it  is  meant  that  two  or 
more  regression  models  are  simultaneously  estimated  under  the  assumption 
that  the  error  terms  are  correlated  within  each  overall  observation. 

When  the  regression  models  are  linear,  then  some  important  facts  are 
known.  If  the  design  matrices  of  the  linear  regression  models  are  all 
equal,  then  the  regression  estimates  are  essentially  the  same  whether  or 
not  they  are  obtained  from  separate  regressions  or  from  one  simultaneous 
regression.  On  the  other  hand,  if  the  intercorrelations  between  the 
predictor  variables  in  the  design  matrices  are  low,  then  the  regression 
coefficients  estimated  simultaneously  have  much  lower  variances  than  the 
regression  coefficients  estimated  from  the  separate  regression  models, 
provided  that  the  correlations  within  each  vector  of  error  terms  are 
high. 

The  usual  impact  acceleration  experiment,  however,  has  the  same  pre¬ 
dictor  variables  for  the  preinjury  as  well  as  for  the  injury  model.  The 
only  reason  that  the  probit  estimates  of  j}j  from  the  "standard”  model  are 
different  from  the  probit  estimates  of  from  the  "modified"  model  is 
that  the  functional  form  of  the  probit  model  is  different  from  the  linear 
preinjury  model.  It  is  suspected  though,  that  considerable  improvement 
in  the  probit  estimate  of  could  be  acheived  if  p  were  large  in  absolute 
value  and  simultaneously  the  predictor  variables  for  the  probit  model 


had  low  correlations  with  the  predictor  variables  for  the  linear 
preinjury  model. 


The  substantiation  of  this  conjecture  would  require  another  Monte 
Carlo  investigation.  Furthermore,  the  feasibility  of  such  an  experimental 
design  needs  to  be  considered.  If  the  conjecture  holds  true,  then  an 
experimental  design  could  be  constructed  to  significantly  improve  the 
probit  estimates  of  .  This  might  be  done  by  first  running  all  test 
subjects  at  acceleration  levels  that  are  low  enough  not  to  cause  injury, 
but  that  are  high  enough  to  give  informative  preinjury  readings. 

"Center"  these  readings  by  subtracting  their  average  values  and  dividing 
by  their  estimated  standard  deviations.  Next  try  to  run  subjects  at 
various  low,  intermediate,  and  high  acceleration  levels  such  that  when 
these  levels  are  "centered"  and  paired,  by  subject,  with  the  low 
preinjury  "centered"  levels  a  pair  of  orthogonal  matrices  results. 

For  example,  if  the  centered  design  matrix  for  the  subjects  run  at 
low,  preinjury  levels  is  given  by  X^  and  the  centered  design  matrix  for 
subjects  run  at  various  low,  intermediate,  and  high  acceleration  levels 
is  given  by  then  X^ 'X^  =  (5.  This  type  of  experimental  design  should 
have  the  property  to  decrease  the  variances  of  the  probit  estimates  for 
every  decrease  in  p.  Furthermore,  moderate  to  high  negative  values  of  p 
should  result  in  substantial  decreases  in  the  MSE  of  the  probit  estimates 
obtained  simultaneously  using  the  preinjury  data  compared  with  those 
obtained  from  the  "standard"  model. 


-9 


IV.  REFERENCES 


1.  Johnston,  J.,  Econometric  Methods.  McGraw-Hill,  New  York,  1963. 

2.  Peterson,  J.  J.,  and  Smith,  D.  E.,  "A  General  Statistical 
Approach  for  Using  Auxiliary  Information  in  the  Development 
of  an  Impact  Acceleration  Injury  Prediction  Model,"  Technical 
Report  No.  112-6,  Desmatics,  Inc.,  August  1980. 

3.  Peterson,  J.  J.,  and  Smith,  D.  E.,  "Computational  Aspects  of 
Incorporating  Auxiliary  Information  into  an  Impact  Acceleration 
Injury  Prediction  Model,"  Technical  Report  No.  112-7,  Desmatics, 
Inc.,  February  1981. 


10' 


UNCLASSIFIED 


security  classification  of  this  page  (Whtn  s««  Eni<»« 


REPORT  DOCUMENTATION  PAGE 


REPORT  NUMBER 


112-9 


GOVT  ACCESSION  NO.I  1 

V'4/OS  */&3\ 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


RECIPIENT'S  CATALOG  NUMBER 


4.  TITLE  (m\d  Submit) 

A  MONTE  CARLO  STUDY  OF  THE  USE  OF 
AUXILIARY  INFORMATION  IN  THE  DEVELOPMENT 
OF  AN  IMPACT  ACCELERATION  INJURY  PREDICTION 
MODEL 


7  AUTHOR  (!) 


TYPE  OF  REPORT  *  PERIOD  COVEREO 


Technical  Report 


PERFORMING  ORG.  REPORT  NUMBER 


contract  or  grant  number^#; 


Dennis  E.  Smith  and  John  J.  Peterson 


N00014-79-C-0128 


9.  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

Desmatics,  Inc. 

P.  0.  Box  618 

State  College,  PA  16801 


<1.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

Office  of  Naval  Research 
Arlington,  VA  22217 


I 


10.  PROGRAM  ELEMENT.  PROJECT,  TASK 
AREA  •  WORK  UNIT  NUMBERS 


NR  207-037 


12.  REPORT  DATE 

September  1981 


IS.  NUMBER  OF  PAGES 

10 


MONITORING  AGENCY  NAME  •  AOORESSf//  dllltttnl  I tom  ContraMnt  Olllct)  IS.  SECURITY  CLASS,  fol  chit  rtporl) 


ISa.  DECLASSIFICATION  'DOWNGRADING 
SCHEDULE 


l«.  DISTRIBUTION  STATEMENT  (ol  Chit  Rtporl) 


Distribution  of  this  report  is  unlimited 


<7.  DISTRIBUTION  STATEMENT  (o I  Iht  mb. Ittcl  tnltrtd  In  Block  30,  II  dllltrtnl  Iron  Rtporl) 


I*.  KEY  WORDS  rConiln 


•<d«  II  nococtory  and  Identity  by  block  numbor) 


Impact  Acceleration  Injury  Prediction 
Auxiliary  Information 
Monte  Carlo  Evaluation 


Z0  ASfT^CT  (Contlnuo  on  rovotoo  »ido  If  ntcatacry  and  identity  by  block  nutnboe) 

This  report  describes  a  small-scale  Monte  Carlo  investigation  of 
procedures  for  incorporating  various  sources  of  auxiliary  information  into 
an  impact  acceleration  Injury  prediction  model.  Parameter  estimates  are 
tabulated  and  compared  for  standard  and  modified  models.  Based  on  the 
results  of  the  investigation,  the  procedures  appear  to  be  helpful  in 
reducing  the  mean  square  error  of  predictions. 

\ 


ssr*  1473 


EDITION  OF  I  NOV  ts  IS  OBSOLETE 


UNCLASSIFIED _ 

SECURITY  CLASSIFICATION  OF  THIS  PAGE  flWi.n  Data  bln  if 


I 


